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Abstract 

We consider the radial wave equation in similarity coordinates within the semigroup formalism. It 
is known that the generator of the semigroup exhibits a continuum of eigenvalues and embedded in 
this continuum there exists a discrete set of eigenvalues with analytic eigenfunctions. Our results show 
that, for sufficiently regular data, the long time behaviour of the solution is governed by the analytic 
eigenfunctions. The same techniques are applied to the linear stability problem for the fundamental 
self-similar solution xt of the wave equation with a focusing power nonlinearity. Analogous to the free 
wave equation, we show that the long time behaviour (in similarity coordinates) of linear perturbations 
around xt is governed by analytic mode solutions. In particular, this yields a rigorous proof for the 
linear stability of xt with the sharp decay rate for the perturbations. 

1 Introduction 

1.1 Motivation 

The focusing semilinear wave equation 

X u -A X = x P (1) 
for x '■ K x R 3 — > R, where p > 1 is an odd integer, exhibits radial self-similar solutions, i.e., solutions of 
the form x(t, x) = (T - t)~ 2/ ^~ 1 ' ) f(\x\/(T - t)) for a function / : R -> R and fixed T > 0. In fact, the 
simplest solution of this type, where / is just a constant, can be obtained by neglecting the Laplacian 
in Eq. ([TJ and solving the resulting ordinary differential equation in t. We refer to this solution as 
the fundamental self-similar solution and denote it by %t- Although self-similar solutions do not have 
finite energy, one may use them together with smooth cut-off functions and finite speed of propagation 
to demonstrate blow up for solutions with smooth compactly supported initial data. This observation 
immediately raises the question how typical such a self-similar blow up is. Does it happen only for 
the very special initial data constructed by the procedure described above or can it be observed for a 
larger set of data? Numerical investigations [T] indicate that the latter is true. Actually, there is a much 
stronger conjecture, namely that the fundamental self-similar solution describes the blow up behaviour 
for generic large initial data. This conjecture is based on numerical investigations for the radial equation. 
In these simulations one observes that the future development of sufficiently large initial data converges 
to the fundamental self-similar solution near the center r = pQ. This indicates that xt has to be stable 
in some sense. We remark that for p = 3 there are also rigorous results in this direction (see [7J, [B], 
.5.). In fact, Merle and Zaag have rigorously proved the full nonlinear stability of a more general family 
of explicit solutions (which includes xt) for the corresponding problem in one space dimension and any 
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p > 1 (see [8], p. 48, Theorem 3). The stability holds in the topology of the energy space. We also 
mention the two recent papers [TO], [pj on interesting consequences of this result. 

In order to analyse linear stability of the fundamental self-similar solution it is convenient to introduce 
similarity coordinates (r, p) defined by r := — log(T— t) and p := s^rj- Since convergence is only expected 
near r = 0, one requires p € (0, 1) which corresponds to the interior of the backward lightcone of the 
spacetime point (t,r) — (T, 0). Transforming Eq. |TJ to similarity coordinates, inserting the ansatz 
X = Xt + <$> and linearizing in <f> yields a rather nasty equation of the form 

1 — o 2 

4>tt + 4>t + 2p4>r P - (1 - p 2 )(f)pp - 2 — - — 4> P - PCo4> = (2) 

where Co > is a constant defined by \t- The first step in a heuristic stability analysis is to look for 
mode solutions, i.e., one inserts the ansatz 4>{t,p) = e Ar u(p). This yields the generalized eigenvalue 
problem 

- (1 - p 2 )u" - 2 1 - P u + 2\pu + [A(l + A) — pco]u = (3) 
P 

which has two singular points at p = and p = 1. A necessary condition for linear stability of \t is the 
nonexistence of mode solutions with ReA > 0. However, it is an entirely nontrivial question what kind 
of solutions of Eq. ((3| one should consider as admissible. In other words, it is not clear what boundary 
conditions one should impose at the singular point p = 1. A basic Frobenius analysis shows that around 
p = 1 there exists an analytic solution and a nonanalytic one where the latter behaves as (1 — p) 1 for 
p —* 1 (we assume noninteger A for simplicity). This shows that the nonanalytic solution becomes more 
and more regular at the backward lightcone as ReA decreases. Hence, if ReA is sufficiently small, there 
is no singular solution which can be excluded a priori. Another difficulty we encounter is the fact that, 
since this is a highly non self-adjoint problem, the nonexistence of unstable modes does not imply linear 
stability. 

The only way to overcome these obstacles is to look for a well-posed initial value formulation for 
Eq. ([2]). It turns out that the machinery provided by semigroup theory can be successfully applied here. 
Very sketchy, one writes Eq. (|2| as a first order system of the form 

jU(r) = LS(t) (4) 

where L is a spatial differential operator which is realized as an unbounded linear operator acting on 
a Banach space. The formal solution of this equation is $(r) = exp(rL)"I>(0) but this does not make 
sense mathematically since L is unbounded. With the help of semigroup theory one is able to construct 
a well-defined one-parameter family S(t) of operators such that the solution of Eq. Q with initial data 
$(0) is given by $(t) = S(t)<&(0). Such a formulation solves the two problems described above. First, 
there exists a well-defined notion of spectrum which implicitly yields the correct boundary condition for 
Eq. ©, and, secondly, one may use abstract results from semigroup theory to obtain growth bounds for 
the solutions. 

1.2 The problem of analytic modes 

For simplicity one may first develop a semigroup formulation for the free wave equation, i.e., Eq. ([2]) with 
Co = 0. This problem has recently been considered [2] and we have shown that there exists a semigroup 
So(t) that yields the time evolution in energy space, i.e., for very rough data. It should be remarked 
that this is an interesting result per se, at least from the mathematical point of view, since the semigroup 
generator is highly non self-adjoint. In fact, it is not even normal and its spectrum has a remarkable 
structure: It consists (essentially) of a continuum of eigenvalues filling a left half-plane in the set of 
complex numbers. We review the corresponding results in Sec. [2] A special subset {0, —I, —2, . . . } of the 
point spectrum consists of eigenvalues with analytic eigenf unctions. From the point of view of semigroup 
theory there is no reason to consider these "analytic eigenvalues" as distinguished. However, in numerical 
evolutions one observes that the asymptotic behaviour (for r — » oo) of solutions is exactly described by 
the analytic eigenvalues and eigenfunctions 0. Therefore, the question is how to explain this behaviour. 
Note that this is not a mere effect of preservation of regularity. In the abstract approach, preservation 

To be precise, this is true only for data that do not have compact support since otherwise Huygens' principle applies. 
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of regularity is expressed by the fact that domains of powers of the generator Lo are invariant under 
the time evolution, i.e., if $(0) G ~D(L k ) for k G N then So(t)$(0) G T>(Lq). But one cannot get rid 
of "nonanalytic eigenvalues" by prescribing data in T>{L^) since any eigenvector of Lo is by definition 
also an eigenvector of Lq. However, in Sec. 13.11 we show that another class of higher Sobolev spaces, 
denoted by Ti 2k , remains invariant under So. A key observation in this respect is a certain commutator 
property exhibited by the generator Lo, see Lemma |3. II below. The spaces Ti 2 are suitable to get rid 
of the continuum eigenvalues and only analytic ones remain. More precise, we show that initial data in 
TC 2k can be expanded in a sum of the first 2k analytic eigenfunctions of Lo plus a remainder whose time 
evolution decays faster than the rest. This result shows in particular that the long time behaviour of 
solutions with smooth initial data is described by the analytic modes as is observed numerically. 

1.3 Application to the semilinear wave equation 

Numerical studies of Eq. ([2} exhibit a very similar behaviour as described above for the free wave equation: 
The large r behaviour of linear perturbations around \t is precisely described by analytic modes, i.e., 
analytic solutions of Eq. The techniques explained above for the free wave equation carry over to 
this problem. We obtain the analogous result (see Theorem 14.11 below) which shows that the long time 
behaviour is indeed given by the analytic modes. In particular, this result yields a rigorous proof for 
the linear stability of the fundamental self-similar solution of Eq. JTJ with the sharp decay rate for the 
perturbation. 

Finally, we remark that many aspects of the problem of analytic modes are related to the work of N. 
Szpak on quasinormal mode expansions for solutions of the wave equation |12j . However, the results in 
[12] have been obtained by very different methods involving the Laplace transform. It is likely that the 
techniques of [12] can also be applied to our problem and this would lead to a very different proof of our 
results. 

1.4 Notations 

To improve readability we write vectors as boldface letters and the components are numbered by lower 
indices, e.g. u = (ui,U2) T . The notation X <— > Y for two normed vector spaces X,Y means that X is 
continuously embedded in Y. When given an inner product on a vector space X we denote the 

induced norm by || • ||x, i-e., || • ||x := \J (-\-)x- The Cartesian product X x Y of two vector spaces X 
and Y with inner products (-\-)x and (-\-)y is implicitly assumed to be equipped with the inner product 
(u\v)xxy '■= (ui\vi)x + (u2\v2)y- For a Banach space X we denote by B(X) the space of bounded 
linear operators on X. For a closed operator L : L>(L) cX->Xwe set Rl(X) := (A — L) _1 whenever 
the right-hand side exists. The resolvent set of L is denoted by p(L) and the point, continuous and 
residual spectra by a p (L), a c (L) and oy(L), respectively (see [2] for the precise definitions). Finally, the 
expression A < B means that there exists a C > such that A < CB. 

2 Semigroup formulation in energy space 

In this section we review results recently obtained by the author [2] on a semigroup formulation of the 
free wave equation in similarity coordinates. We define similarity coordinates (r, p) as explained in the 
introduction by r := — log(T — t),p:= tj?—^ and consider the radial wave equation on (3 + 1) Minkowski 
space, 

2 7 

4>tt — Iprr T/V = 0. 

r 

Substituting ip(t,r) := rip(t,r) yields 

Iptt ~ Iprr = 

with the boundary condition ^(£,0) = for all t. We write this equation as a first order system 
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where $ := (ipt,ip r ) T . Changing to similarity coordinates we obtain 

dr* = ( 7 \ ) d P * (5) 

where $(r, p) := *(T - e~ T , pe~ r ). 

Let ft := L 2 {0, 1) x L 2 (0, 1), V(L ) := {u G C^O, 1] x C* 1 ^, 1] : ui(0) = 0} and 

T „( n \.-( -P u 'i(P) + U '-2(P) \ 
L ° U(P) -1 <{p)-pu' 2 {p) )■ 

I/O : T>(Lo) C ft — > ft is a densely defined linear operator on the Hilbert space ft. An operator formulation 
of Eq. (0 is given by 

4-*(f) = L <fr(r) 

for a strongly differentiable function $ : [0, oo) — > ft. We have the following result [2J. 
Theorem 2.1. TVie operator Lq is closable and its closure Lq generates a strongly continuous one- 
parameter semigroup So : [0, oo) — > B(ft) satisfying ||S'o(t)||b(?<) < e5T / or Q ^ r > 0. 

TTie spectrum of Lo is given by a p (Lo) = {A G C : ReA < -}, cr c (Lo) = {A G C : ReA = 
cr r (L ) = 0- 

3 Semigroup formulation for more regular data 
3.1 Invariance of higher Sobolev spaces 

We show that a certain class of higher Sobolev spaces is invariant under the semigroup So- For k G No 
we set 

ft 2fc := {u G H 2k (0, 1) x H 2k (0, 1) : 4 23) (0) = u (2]+1) (0) = 0, j G No, j < fc} 

and define an operator D 2 : ft 2 — > ft by D 2 u := u". We have ft = ft and equip ft 2fe with the inner 
product (u|v) H 2k := (u|v)-h + (D 2k \\\D 2k v)-n. The following lemma summarizes elementary properties. 

Lemma 3.1. 1. TL 2k is a Hilbert space. 

2. Jl 2( - k+r > is a dense subspace ofTt 2k and the inclusion 7-^ 2 ( fe+1 ' c Ti 2k is continuous. 

3. The operator D 2 satisfies D 2 H 2{k+1) C H 2k . 

4. We have H 2(k+1) C T>(L ) and L ft 2(fe+1) C H 2k . 

5. D 2 and Lo satisfy the commutator relation D 2 Lou = Lo-D 2 u — 2D 2 u for all u G ft 4 . 

Proof. The proof is straightforward by inserting the definitions and using well known properties of 
Sobolev spaces. □ 

As usual we define the part L , k of Lo in H 2k by V{L ,k) := {u G V(L ) n Ti 2k : Lou G H 2k } and 
Lo,fcU := Lou. We show that Lo,fc generates a semigroup on ft 2fc . 

Proposition 3.1. The operator Lo,k generates a strongly continuous one-parameter semigroup So } k : 
[0,oo) -> B(H 2k ) satisfying \\So,k\\ B (n™) < e ' T ■ 

Proof. By Lemma 13.11 we immediately observe that Lo,fc is densely defined since ft 2 ^ 1 ' c V(Lo,k). 

Let (uj) C 2?(L ,fc) with Uj — > u and L ,fcUj — » f both in H 2k . Since H 2h ft (Lemma [3TTJ 
this implies ir, — » u, Lou., — > f in ft and by the closedness of Lo we conclude u G T>(Lo) n H 2k and 
Lou = f G ft 2fe which shows u G T>{Lo,k) and we have proved that Lo.fc is closed. 

By using the commutator relation from Lemma |3. II and integration by parts (cf. |2j) we obtain 

Re(L ,fcu|u) W 2fc = Re (^(L u\u) n + (L D 2k u\D 2k u) n - 2fc||D 2fc u|| 2 ^ < -||u||^ 2 «r 
for all u G fi 2 ^ k+1 > and by a density argument this estimate holds in fact for all u G T)(Lo,k)- 
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Let f G 7i 2k n C°°(0, l) 2 and define F(p) := Mp) + pf 2 (p) + /> u 2 (p) := ^ f p F(£)d£ 

and Mi(p) := pu 2 (p) — f P f2((,)d£ i . Then the Taylor series expansion for u\ around p = up to order 
2k — 1 contains only odd powers of p whereas the analogous series for u 2 up to order 2k contains only 
even powers of p. This shows that u satisfies the appropriate boundary conditions at p = and we 
conclude that u G TC 2k (1 T>{Lq). Furthermore, a direct computation yields (1 — Lo)u = f which shows 
that u G T>(Lo t k) and 1 — Lo,fe has dense range. 

Invoking the Lumer-Phillips Theorem (see e.g. [3], p. 56, Theorem 4.2.6) finishes the proof. □ 

Based on this result we are able to conclude the invariance of 7i 2k under the semigroup So. 

Lemma 3.2. The space H 2 is Lo -admissible, i.e., it is an invariant subspace of So(t), r > 0, and 
the restriction of So(t) to TC 2k is a strongly continuous semigroup on H 2k satisfying ||So(t)u|| W 2)< < 
e5 T ||u|| W 2fc for all u G H 2k and r > 0. 

Proof. Let f G H 2k and A G p(Lo). Proposition 13.11 implies that there exists a u G X>(Lo,fc) such that 
(A — Lo,fc)u = f. However, since Lo,fc C Lo, we have (A — Lo)u = f and thus, 7?L (A)f = u G Ti. . This 
shows that RL {X)TC 2k C 7i 2k . By Lemma |3. II the embedding fi 2k C TL is continuous and therefore, the 
claim follows from Proposition 13.11 and the theorem on admissible spaces (see e.g. [11] , p. 123, Theorem 
5.5). □ 



3.2 Decomposition 

We improve the growth estimate HSoM \ n 2h \\B(n 2k ) ^ e ^ T by a decomposition of the initial data space 
Ti 2k . Let Af denote the set of all u G TL 2k such that D 2k u = 0. Af is a finite dimensional (and hence 
closed) subspace of H 2k ■ Thus, there exists a bounded orthogonal projection P G B(Tt 2k ) associated to 
it. We immediately obtain ||5'o(r)u|| H 2fc < ||S'o(r)Pu|| H 2fc + \\S a (r)(I - P)u\\ H 2k for any u G Tt 2k . In 
order to keep things simple, we analyse the two parts j|S'o(T)Pu|| 7i 2): and \\Sq(t)(I — P)u\\ n 2k separately. 

Lemma 3.3. The subspace Af is spanned by 2k analytic functions u(-, Xj), j = 1,2, ... ,2k where each 
u(-, Xj) is an eigenfunction of Lo with eigenvalue Xj — —j + 1. 

Proof. For j = 1,2, ... ,2k set u(p, Xj) := (l-p) 1_Aj -(l+p) 1 ~ Xj , ui(p, Xj) := pu' (p, Xj) + (Xj — l)u(p, Xj) 
and u 2 (p, Xj) := u'(p, Xj). Then u(-, Xj) is an eigenfunction of Lo with eigenvalue Xj as a straightforward 
computation shows. Observe that u(-,Xj) is an odd function (binomial theorem, all even powers cancel) 
and therefore, ui (■, Aj) is odd and u 2 (-, Xj) is even. This shows that u(-, Xj) G Ti 2k for each j. Note further 
that Ui(-,Xj) and u 2 (-,Xj) are polynomials of degree strictly smaller than 2k and thus, D 2k u(-, Xj) — 
for each j and we conclude u(-, Aj) G Af. A function in TL 2k satisfies 2k boundary conditions and thus, the 
space Af, which consists of u G Ti 2k with D 2k u — 0, is 2fc-dimensional. However, the 2k eigenfunctions 
u ('> 3 ' = 1> 2, • . • , 2k belong to Af and they are linearly independent. □ 

Lemma [3. 31 shows that we can calculate the time evolution So(t)Pu explicitly: Since Pu £ Af we can 
expand it in terms of analytic eigenfunctions u(-, Xj) of Lo, 

2k 

Pu = ]Pc 3 -u(-,Aj), 

3=1 

where Aj = — j + 1 and c 1 ,c 2 , . . . , c 2 k G C are the expansion coefficients. Hence, we obtain 

2k 

5o(r)Pu = ^ C ,^ T u(-,A,). 

3=1 

In order to estimate Sq(t)(I — P)u we need a few preparations. First, we state the following general 
observation. 

Lemma 3.4. Let X,Y be Hilbert spaces with X C Y and A : X — > Y a linear bounded surjective 
mapping. Then there exists a constant c > such that \\Au\\y > c||u||js: for all u G {ker A)^ . 
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Proof. By the boundedness of A we conclude that ker A is closed and thus, we have the decomposition 
X = ker A © (kerA) 1 . Hence, the mapping ^(kcrA)^ i s bijective. The closed graph theorem implies 
that its inverse is bounded as well and thus, there exists a c > such that > c||u||x for all 

uG(ker^) ± . □ 

Now we are able to prove the desired estimate for So(t)(I — P)u. 
Proposition 3.2. 1. The subspace Af C H 2k is invariant under the semigroup So- 

2. The mapping uh ||D 2fe u||H defines a norm on Af ± which is equivalent to \\ ■ W^k. 

3. For u G Af 1 " we have the estimate ||5o(t)u[| W 2* < e { ^~ 2k)T ||u|| W 2* for all r > 0. 

Proof. 1. Let f G Af ± (orthogonal complement in H 2k ), f 7^ 0, A G p(Lo) and set u := Rl {X){. Then 
we have u G Ti. 2k (cf. proof of Lemma |3,2[) . Suppose u ^ Af ± , i.e., u G Af. Then the commutator 
relation from Lemma I3TT1 yields D 2k i = D 2k (\ - L )u = L D 2k u = which implies f G TV and 
we infer f = 0, a contradiction. Thus, we have R.L (\)Af ± C Af ± . Since Af ± is closed in H 2 and 
So{T)\ n 2k defines a semigroup on TC 2k by Lemma 13.21 the claim follows from [TT], p. 121, Theorem 
5.1. 

2. The mapping D 2k : TL 2k — > Ti is linear, bounded and surjective. According to Lemma 13.41 we have 

< l|0 2fc u|| H for all u £ Af 1 - = (kerD 2k ) ± . Trivially, we have ||£> 2fe u|| H < \\u\\ n 2k. 

3. We equip Af 1 with the inner product (u|v)jvx := {D 2k u\D 2k v) n - Then Af x is a Hilbert space 
and So(t)\j^± defines a semigroup on Af 1 by assertions 1 and 2 from above. The generator of this 
semigroup is L , the part of Lq in Af 1 " , and it satisfies 

Re(L 0M ± U \uy± =Re(L D 2k u\D 2k u)u-2k\\u\\ 2 ul _ < Q - 2fc) ||u||^ 

for u £ T>(L j^±) where we have used the commutator relation from Lemma |3. II iterativelv. This 

estimate implies ||S'o(r)u|| A /-j. < e < -3 _2fc ' T ||u|| jV - i . However, by assertion 2 above, the norm || ■ \\j^± 
is equivalent to || ■ \\ n 2k and we arrive at the claim. 

□ 

Proposition 13.21 implies that \\Sq(t)(I - P)u\\ n2 k < e^~ 2k)r ||u|| w aii since (I - P)u G N 1 and this 
completes the investigation of the free wave equation in similarity coordinates. We end up with the result 
that the time evolution can be estimated as 

2k 

l|So(r)u|| wafe <J2^ jT + ^' 2k)T Mn- 
3=1 

for initial data u G H 2k and Xj = + 1 if we assume the eigenfunctions u(-,Aj) to be normalized. 
This completely answers the question of the role of the analytic modes. If the initial data are sufficiently 
regular then the long time behaviour of the solution is dominated by the first analytic eigenmodes. This 
is exactly what is observed numerically. We also emphasize that this result implies a certain completeness 
property of the analytic modes. Sufficiently regular initial data can be expanded in a sum of analytic 
modes plus a remainder which decays faster. We summarize the results in a theorem. 

Theorem 3.1. Let u G Tt 2 . Then there exist constants ci, . . . , C2k G C and a function g G TL 2k such 
that 

2k 
3=1 

and ||S'o(r)g|| W 2fc < e'2 -2fc ' T ||g|| H 2fc where u(-,Xj) are normalized analytic eigenfunctions of Lq with 
eigenvalues Xj = — j + 1. In particular, we have 

2k 

HS>(t)u||„» < E + ^- 2k)T \\u\\ n 2k 

3=1 

for all t > 0. 
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4 Application to the semilinear wave equation 



We apply the previously obtained results to the linear stability problem for the fundamental self-similar 
solution of Eq. fTJ. To this end we construct a semigroup S acting on Ti. 2k that describes the time 
evolution of linear perturbations of the fundamental self-similar solution. Througout this section we 
restrict ourselves to k £ N since the case k = has already been investigated in [2J. 



4.1 Operator formulation 

Let V £ B(H) be defined by 

L'u(p) := 







where the constant Co is given by the fundamental self-similar solution \t- Note that L' leaves H 2k 
invariant as an odd-even argument easily shows. Furthermore, D 2k and L' commute, i.e., D 2k L'u = 
L'D 2k u for all ueH 2k . 

An operator formulation for the linear stability problem Eq. © is given by 

jU(r) = L*(t) (6) 

where L := Lo + L' and the Bounded Perturbation Theorem (see e.g. [3]) immediately yields the existence 
of a semigroup S : [0, oo) -> B(H) satisfying ||5(r)[[ B (w) < e ( s +peo)T for all r > and the solution $ of 
Eq. © is given by $(r) = 5(t)$(0). 

4.2 Invariant subspaces 

Lemma 4.1. TTie space Tl 2k is L-admissible, i.e., it is an invariant subspace of S(t), t > 0, and 
the restriction of Sir) to TL 2k is a strongly continuous semigroup on H 2k satisfying \\S(r)u\\ H 2k < 
e (i+pco)T|| u || w2fc j or all u £ n 2k and T>0 

Proof. The part of L in H 2k is given by L , k + L'\ n2k since L'H 2k C H 2k . Thus, Proposition 137X1 and 
the Bounded Perturbation Theorem show that the part of L in TL 2k generates a strongly continuous 
one— parameter semigroup on Tt 2k with the same growth bound as 5*. Applying the same argument as in 
the proof of Lemma [372] yields the claim. □ 



As before, we denote by M the space of all u £ TL 2k such that D 2k u = and by M the orthogonal 
complement of Af in H 2k . 

Proposition 4.1. The subspace Af ± is invariant under the semigroup S. Furthermore, we have the 
estimate 

I|S(t)u||„„ <e^+™~ 2k ^\\u\\ n2k 

for all u £ Af ± and r > 0. 

Proof. Let A > | + pc . Then R L (\)H 2k C H 2k since H 2k is ^-admissible by Lemma ITTI (cf. [TT] , p. 
123, Theorem 5.5). Set u := R L (X)f for an f £ A/' ± , f / 0, and suppose u g A/' ± . Then u £ jV and 
the commutator relations show D 2k f = D 2k {\ — L)u = which implies f £ M and we infer f = 0, a 
contradiction. Hence, Rl{\)M^ C Af ± and, since Af ± is closed in Tt 2k , the claim follows from [llj . p. 
121, Theorem 5.1. 

We define an inner product on M 1 - by Hv)^ := {D 2k u\D 2k v) n . The induced norm || ■ is 
equivalent to || • \\ H 2k on A/" x by Proposition 13.21 and thus, JV ± equipped with (-\-)j^± is a Hilbert space 
since J\f ± is closed in TL 2k . The restriction S(t)\j^± defines a semigroup on J\f ± and its generator is the 
part of L in Af ± , denoted by Lj^± . The generator satisfies 

Re ((L C 2fe u|.D 2fe u)H + {L'D 2k u\D 2k u) n ) - 2k{D 2k u\D 2k u) n 

< Q+pco-2*) ||u||V 
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for all u £ V(L A f±) where we have used the commutator relation from Lemma 13.11 and the fact that 
V and D 2k commute, as already remarked. This, however, implies the growth estimate ||S , (r)u|| jV -j. < 
e ( 2+ pc o- 2fc ) T ||u|| jV - i for all u £ A/' ± and r > and by the equivalence of the norms || ■ \\ n 2k and || ■ \\ M ± 
(Proposition 13. 2( l we arrive at the claim. □ 

Note that the growth estimate given in Proposition 14. II is certainly not optimal, however, we do not 
make any attempts to improve it since we can make k arbitrarily large. 



4.3 Analytic modes 

The point spectrum of L can be calculated by solving a second order ODE. To this end we recall the 
definition of the operator T(A) from [2]. The domain of T(A) is D(T(A)) := {u £ H 1 ^,!) : u £ 
H? oc (0, l),t(X)u £ L 2 (0, l),u(0) = 0} and T(A)u := t(X)u where 

t(X)u(p) := -(1 - p 2 )u"(p) + 2Xpu'(p) + [A(A - 1) - pc ]u(p). 

Then A £ o~ p (L) if and only if dimkerT(A) = 1 (see 2 , Proposition 2). 

Observe that t(X)u = corresponds exactly to Eq. ([3]) in view of the substitution u(p) i— > pu(p). 
Thus, as explained in the introduction, the semigroup approach implicitly yields the correct boundary 
condition for the generalized eigenvalue problem Eq. @ that defines mode solutions. 

The general solution of t(X)u = can be given in terms of Legendre functions and hence, the point 
spectrum can be calculated explicitly. In [I] and H] it has been shown that t(X)u — has a nontrivial 
analytic solution if and only if A = A* where A+ := 1 + — 2j and XJ := — -4rj ~ 2J for a j £ No- We 
will refer to A^ as analytic eigenvalues. Moreover, it follows from [I] that the analytic functions it(-, A*) 
satisfying t(Xf)u(-, A*) = are in fact odd0 polynomials of degree 2j + 1. 

Very similar to the free wave equation, the subspace TV is again spanned by analytic eigenfunctions 
of L. 

Lemma 4.2. The subspace M is spanned by 2k analytic functions u(-, A*) for j = 0, 1, 2, . . . , k — 1 where 
each u(-, A*) is an eigenfunction of L with eigenvalue A* . A^ = 1 + — 2j and XJ = — ^rr — 2j. 

Proof. The space M is 2fc-dimensional as already remarked in the proof of Lemma 13.31 Let u(-, A*) 7^ 
satisfy t(Xf)u(-, Xf) = for j = 0, 1, . . . , k - 1 and define u(-, Xf) by m(p, Xf) := pu'(p, Xf) + (Xf - 
l)it(p, A^) and «2(p, A*) := u'(p, Xf). Then u(-, A*) £ TL 2k n TJ(L) since ui(-, A*) is an odd polynomial 
and U2(-,Xf) is an even polynomial. A straightforward computation shows Lu(-,Xf) = Xfu(-,Xf) 
and thus, u(-,A*) is an eigenfunction of L with eigenvalue A*. The functions ui(-,Xf) and U2(-,Xf) 
are polynomials of degree strictly smaller than 2k and therefore, D 2fc u(-,A*) = which shows that 
u(-, Xf) £ M for all 3 = 0, 1, 2, . . . , k - 1. However, since the u(-, Xf) are eigenfunctions with different 
eigenvalues they are linearly independent and they form a set of 2k linearly independent functions in the 
2fc-dimensional space TV. □ 



4.4 The time evolution of linear perturbations 

We have collected all the necessary preliminaries to conclude the analogous result to Theorem 13.11 
Theorem 4.1. Let u £ H 2k . Then there exist 2k constants Cq , ... ,c^_ 1 £ C and a function g £ Ti 2k 
such that 

fc-i 

u = E ( C M-' x t) + c 7 u (-> A 7)) + s 

j=0 

and ||S'(r)g|| H 2fc < e^ +vc °~ 2k ^ T \\g\\ n 2k where u(-,Aj = ) are normalized analytic eigenfunctions of L with 
eigenvalues XJ = 1 + — 2j and XJ — — — 2j. In particular, we have 

\\S(r)u\\ H2k <E(c+e A + T + C -e A ^) + e^"- 2 ^ \\u\\ n2k 
3=0 



2 Note that in our convention the whole spectrum is shifted to the right by — compared to [T]. 

3 In [l] the corresponding analytic "eigenfunctions" are even polynomials of degree 2j but, according to our convention, one 



has to multiply them by p. 
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for all t > 0. 

By making k sufficiently large we infer that the long time behaviour of smooth perturbations is 
governed by the analytic modes and this is exactly what is observed numerically. 

Furthermore, the largest eigenvalue Ay = 1 + is known to emerge from the time translation 
symmetry of the original problem (cf. [2]). This apparent instability is merely an effect of the similarity 
coordinates and it is therefore called the gauge instability. Thus, for studying the question of linear 
stability we only allow perturbations with Cq — (notation as in Theorem 14. i.e., perturbations such 
that the gauge instability is not present. Theorem |4T] shows that the time evolution of sufficiently regular 

perturbations u 6 7i 2k with Cg — decays as ||5 , (r)u|| H 2fc < e p- 1 T \\u\\ n 2k for r — > oo and this estimate 
is clearly sharp. Hence, the decay is exactly described by the largest analytic eigenvalue apart from the 
gauge instability. We conclude that the fundamental self-similar solution for the wave equation with a 
focusing power nonlinearity is linearly stable. 
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